Enhanced Pairing in the "Checkerboard" Hubbard Ladder 
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We study signatures of superconductivity in a 2-leg "checkerboard" Hubbard ladder model, de- 
fined as a one-dimensional (period 2) array of square plaquettes with an intra-plaquette hopping t 
and inter-plaquette hopping t' , using the density matrix renormalization group method. The highest 
pairing scale (characterized by the spin gap or the pair binding energy, extrapolated to the thermo- 
dynamic limit) is found for doping levels close to half filling, U ~ 6t and t /t ~ 0.6. Other forms 
of modulated hopping parameters, with periods of either 1 or 3 lattice constants, are also found to 
enhance pairing relative to the uniform two-leg ladder, although to a lesser degree. A calculation 
of the phase stiffness of the ladder reveals that in the regime with the strongest pairing, the energy 
scale associated with phase ordering is comparable to the pairing scale. 
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I. INTRODUCTION 

The much debated theoretical issues related to the 
"mechanism" (i.e. microscopic origin) of high tempera- 
ture superconductivity are often ill-defined. One related 
question to which unambiguous answers are possible is: 
For a given class of models, what values of the param- 
eters are optimal for superconductivity? Of course, if 
one can make predictions about models, the same in- 
sights might provide guidance in the search for materials 
with improved superconducting properties. Two specific 
questions we would like to address are: 1) In the case in 
which superconductivity arises directly from the repul- 
sive interactions between electrons, how strong (in units 
of the bandwidth) are the optimal interactions for super- 
conductivity? 2) Is there an "optimal inhomogeneity" for 
superconductivity 1 , in the sense of a complex (but still 
periodic) electronic structure with multiple orbitals per 
unit cell? An obvious difficulty with this program is that, 
in most cases, we do not know how to compute the tran- 
sition temperature of the relevant models in a controlled 
manner, so as to test the predictions of theory. 

In this context, we use density matrix renormalization 
group (DMRG)£ to numerically compute the supercon- 
ducting correlations of the two-leg Hubbard ladder (ex- 
trapolated to infinite length) as a function of the strength 
of the Hubbard interaction, U, and for various periodic 
patterns of the hopping matrix elements. The ID charac- 
ter of the system studied is what permits us to obtain an 
accurate solution of this problem. However, the same ID 
character implies that no non-zero critical temperature 
is possible, so in assessing the optimal conditions for su- 
perconductivity, we are forced to use other energy scales 
in the problem, especially the spin-gap, AE S , the pair- 
binding energy, AE p , and the superfluid hclicity modu- 
lus, p c . We find that: 1) The optimal value of U is gen- 
erally U ~ 6t where 6t is the total bandwidth of the uni- 
form ladder. This result agrees with previous studies 3 -"^ 
of various ladder systems. It is also consistent with in- 
ferences made on the basis of exact diagonalization^ and 



dynamical cluster quantum Monte- Carlo^ studies of the 
2D Hubbard model, where U ~ 8t (i.e. the 2D band- 
width) was found to be optimal. 2) For the checkerboard 
pattern with 4 sites per unit cell shown in Fig. lb, the 
optimal conditions occur for an intermediate degree of 
inhomogeneity, t' jt ~ 0.6 — 0.7, where t is the hopping 
matrix within a square and t' is the hopping matrix be- 
tween squares. This tends to corroborate inferences made 
previously on the basis of exact diagonalization studies 5 
of the 2D "checkerboard-Hubbard model." 3) A quali- 
tatively similar enhancement of superconductivity is ob- 
served for the other periodic versions of the model with 
2 or 6 sites per unit cell shown in Figs, la and lc, re- 
spectively, although in these cases the magnitude of the 
effect is smaller and the optimal condition occurs with 
values oit'/t closer to 1. 

The observation that certain patterns of spatial sym- 
metry breaking can coexist with superconductivity (or 
even strongly enhance it), while others do not, is also 
reminiscent of recent results obtained using DMRGi and 
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FIG. 1: Schematic representation of the "inhomogeneous" 
Hubbard ladders considered in the present paper: a) The pe- 
riod one "dimer" ladder; b) the period two "checkerboard" 
ladder; c) the period three ladder. As discussed below Eq. 
the solid and dashed lines represent, respectively, hopping 
matrix elements t and t' . 
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the dynamic cluster approximation 8 . In the first of these 
calculations, the inhomogeneity (in the form of stripes) 
occurs spontaneously, while in the second it is imposed 
externally. As we were completing this work, a contractor 
renormalization (CORE) study of the checkerboard Hub- 
bard model in a 2D geometry was presented in Ref. |9|, 
extending earlier CORE results for the uniform 2D Hub- 
bard modeliS. Finite size effects were found to be large 
for t' > 0.8i, but in the smaller t' regime, where these 
effects are relatively small, the results of this new study 
are completely consistent with those of the earlier exact 
diagonalization studies^, and lead to conclusions concern- 
ing the optimal conditions for superconductivity that are 
similar to those obtained in the present ladder study. The 
CORE method was also used to study ladders, albeit con- 
siderably shorter than those studied in the present paper, 
and again the results obtained are fully consistent with 
the present results. 



II. THE MODEL 

We consider the repulsive U Hubbard model defined 
on a (spatially modulated) two-leg ladder 



*XMi,<r c J.2, CT + h.c.) + U X/'.-a.-'VA.. (!) 
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Here cj A g creates an electron on rung j ~ 1 , . . 
of chain A = 1, 2 with spin polarization a — ±, L is 
the length of the ladder, U > is the repulsion between 
two electrons on the same site, the density operator is 

%-,A,<7 = c] Xa c hX M, and n = (2L)" 1 Ea,>i,A, ff ) is 
the mean number of electrons per site. The much stud- 
ied homogeneous Hubbard ladder corresponds to the case 
tjj+i — t' for all j, although it is worth noting that for 
t' <C t, this model can also be viewed as a coupled array 
of Hubbard-dimers. The "dimer ladder" is shown in Fig. 
QJ» The "checkerboard ladder" in Fig. [T]d has i2j,2j+i = t 
and t2j+i.2j+2 = f < t. The "period three" ladder in Fig. 
QJ has tsj^j+i = £3^+1,3^+2 = t and i3j+2,3j+3 = t' <t. 

The thermodynamic limit is accessed by computing 
quantities for various lengths, and then using finite size 
scaling analysis to extrapolate to 1/L — > 0. 
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III. EFFECTIVE FIELD-THEORY 

The uniform two-leg Hubbard ladder with n ^ 1 but 
still not too far from n = 1, is well known, on the ba- 
sis of weak coupling RG 11 , bosonzationi^, and DMRG 3 
approaches, to be in a Luther-Emery phase character- 
ized at low energies by a spin-gap, AE S (defined in Eq. 
HI below) and a single, gapless acoustic "charge" mode 
which propagates with speed v Cl and whose long-range 



(power-law) correlations are determined by a single Lut- 
tinger parameter, K c . The Luther-Emery liquid can be 
thought of as a ID version of a superconducting state in 
the sense that it has a non-vanishing supernuid stiffness 
(see Eq. below), and, for K c > 1/2 and T < AE S , it 
has a divergent superconducting susceptibility, 



X ~ Xo 



(2-l/Kc) 



(2) 



where v c is the charge velocity and a is a lattice constant. 
In the single chain realization of a Luthcr-Emcry liquid, 



Xo 



aAE, 



(3) 



where £ s — v s /AE s is the spin-correlation length and v s 
is the spin-velocity. For a multicomponent system, the 
corresponding expression for \q is somewhat more com- 
plicated, as there may be multiple scales (e.g. multiple 
spin-gaps) associated with the gapped modes. However, 
Xo remains a monotonic, approximately linearly increas- 
ing function of AE S . 

Perhaps not surprisingly, we will see that the inho- 
mogeneous Hubbard ladders we have studied are also 
Luther-Emery liquids with K c > 1/2. Thus, in address- 
ing the "mechanism of superconductivity," the primary 
purpose of our DMRG calculations is to determine the 
dependence of v c , K c , AE S and £ s on microscopic pa- 
rameters. 

The pair binding energy AE p corresponds to creating 
two spatially-separated spin-1/2 quasiparticles. Since the 
spins for these quasiparticles can either add to S = 
or 1, we must have AE S < AE p . If the residual inter- 
actions between quasiparticles are repulsive, we expect 
AE S = AEp. Conversely, if the interactions between 
quasiparticles are attractive, a neutral spin-1 "exciton" 
is formed, which has lower energy than two far-separated 
quasi-particles, and hence AE S < AE p . The latter be- 
havior has been found previously in DMRG calculations 
on the uniform t — J ladder—. 



IV. ENERGY SCALES 

The spin-gap, AE S , is the difference between the 
ground-state energies of the system with spin S — 1 and 
S = 0: 

AE S =£ (S=1,2N)-S (S = 0,2N), (4) 

where £o(S,N) is the spin S ground-state energy of the 
N electron system. 

Similarly the pair-binding energy, AE p , is defined as 

AE p = 2£ (^,2N + 1) - £ (0, 2N) - £ (0,2JV + 2). (5) 

Were we computing these quantities in a BCS supercon- 
ductor, then in the thermoydnamic limit, both these en- 
ergies would be equal to twice the minimum gap A m j n 
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lim AE S = lim AE p = 2A min . 

Thus, it is intuitively reasonable to associate these energy 
scales with a mean-field estimate of the superconducting 
critical temperature, T C MF ~ AE s /4. Of course, since 
the ladder is a fD system, the actual T c = 0. 

While it may be reasonable to interpret AE S and/or 
AE p as measures of a pairing scale in the problem, in 
order to address the growth of superconducting correla- 
tions it is ultimately necessary to consider the helicity 
modulus, which governs the energetics of superconduct- 
ing phase fluctuations: 



Pc 



v c K c 
2ir 



= lim 



\r d2£ ° 






d4> 2 


0=0 



(6) 



where, in this case, the ground-state energy is computed 
in the presence of pair-fields applied to the two ends of 
the system with a relative phase twist (p. 

In 2D, the relative importance of phase and pair- 
breaking fluctuations can be assessed^ by considering 
the ratio of the phase stiffness (which has units of en- 
ergy) to the pairing gap. However, in ID, p c has units 
of a velocity, so defining an energy scale, AEg, charac- 
teristic of the phase fluctuations requires introducing a 
length scale in the problem. The only emergent length 
scale is £ s , in terms of which we define 



AE e = ir Pc /Zs = R AE S 



(7) 



Here R = AEg/AE s is the dimensionless ratio of the 
phase ordering and pairing scales. 

To appreciate the significance of this ratio, consider 
its value for the attractive Hubbard chain in various lim- 
its. The ID version of a BCS limit, in which there is 
a single characteristic energy/temperature scale, A s ~ 
exp[— 7ru s /a|?7|], is realized in the limit \U\ -C 1 where, 
up to corrections of order U/t, v s = v c and K c — 1, so 
R = v c K c /2v s — l/2+0(U/t), i.e. both mesoscale phase 
coherence and pairing correlations onset at a tempera- 
ture of the order of T pa j r ~ AE S /A. Conversely, R — >• 
as \U\/t — > oo; for large U, a spin pseudo-gap opens when 
T ~ Tp a i r = \U\/2, with a second crossover from a largely 
incoherent paired state to a coherent Luther-Emery liq- 
uid occurring at a temperature Tg ~ AEg oc t 2 /\U\, well 
below Tp a i r ■ A similar dichotomy exists in the two- leg 
repulsive U Hubbard ladder, where R — > as the doped 
hole concentration, x — > 0, while R ~ 1 at larger val- 
ues of x where the spin-gap is significantly suppressed 
relative to its value at x = 0. In the small x case, the 
doped holes can be thought of as a dilute gas of charge 
2e bosons at temperatures small compared to T pa j r , but 
the phase coherence scale is much smaller and vanishes 
as x — > 0. 

With these examples in mind, we identify the case 
R ~ 1 with the ID version of the "BCS-like limit" in 
which there is a single crossover temperature Tp a ; r which 



separates the "normal" (multicomponent Luttinger liq- 
uid) high temperature regime from the low temperature 
regime in which substantial mesoscale superconducting 
order has developed. Conversely, if R -C 1, two dis- 
tinct crossover scales characterize the evolution from the 
normal state: a first, high temperature crossover, T parr , 
characterized by the opening of a spin pseudo-gap, and a 
lower crossover temperature, Tg ~ AEg /4, which can be 
viewed as the scale at which the liquid of bosonic pairs 
begin to exhibit substantial local phase coherence. 

The most direct and efficient way to compute £ s from 
DMRG is to apply a staggered Zeeman field to one end 
of the ladder, j = (thus locally breaking spin-rotational 
symmetry) and then measure the decay of the magneti- 
zation as a function of distance down the ladder. In a 
spin-gapped phase, we expect 

a 

- cos[Qj + 0q] exp[-|j|a/£ s ]. (8) 

In the limit of an asymptotically small spin-gap, Q — 
2kp, but for larger gaps it may depend not only on n but 
on U/t as well. To be explicit, we therefore define the 
spin correlation length as 



Ej \J M(j)\ 

E,|m(j)| ■ 



(9) 



It turns out that Eq. [5] is relatively difficult to imple- 
ment to obtain quantitatively reliable results for p c using 
DMRG. However, it is possible- to compute p e by sepa- 
rately calculating v c and v c /K c from quantities that are 
more straightforwardly computed using DMRG. From 
the bosonized field theory, we can identify the inverse 
compressibility of the ladder with the ratio . In turn, 
in all circumstances relevant to the present calculation^, 
the compressibility is related to the energy to add or re- 
move a singlet pair of electrons from the ladder: 



— = lim L 

K L—>oo 



So (0, 2N + 2) + £ (0, 2N - 2) - 2£ Q (0, 2N) 

4 ' 
(10) 

An independent measurement of v c can be obtained by 
calculating also the energy of the first excited state, 
£ i (S, N) according to 



v c = lim - [£ x (0, 2N) - £ a (0, 2N)} 
We then compute the helicity modulus as 

V 2 K 



Note that this procedure also gives us 



K r = —KV r . 

2 



(11) 



(12) 



(13) 



An alternative way to obtain K c is by measuring the 
amplitude of the "Friedel-like" density oscillations which 



exhibit a power-law decay as a function of distance from 
the edge of the system. For long systems, the density 
near the center of a length L ladder takes the forrn^ 



cos[27rn(j - L/2)} 



(14) 



Therefore, by measuring the amplitude of the density 
oscillations Aqbw vs. L and plotting log(-AcDw) vs. 
log(L), K c can be obtained. Whenever possible, we have 
calculated K c using both Eq. [13] and Eq. [Ml and found 
that the two values agree with each other to within about 
10%. 



V. DMRG RESULTS 

We have computed ground state properties for lad- 
der systems for various values of n, t' /t, and U jt using 
DMRG. We have kept up to m — 2400 states and ex- 
trapolated our results to zero truncation errors. As is 
well known, 1 - ground state energies (as well as one-point 
correlation functions — ) can be extracted with great accu- 
racy in this way. Results have been obtained for system 
sizes from 2 x 16 up to 2 x 64, and then extrapolated to 
the thermodynamic limit (1/L — > 0) using a finite size 
scaling analysis. From the extrapolated values, we have 
extracted AE S , £ s , AE p , p c , and K c , as described above. 

In Fig[2l we show AE S for fixed U/t = 8 as a function 
of t'/t for n = 0.9375, 0.875, and 0.75. It has previ- 
ously been shown that, for the uniform Hubbard two- 
leg ladder with n = 0.875, the spin-gap is maximal for 
U = 8t£ Note that the value of AE S rises from its value 
for the uniform ladder as t'/t is reduced below t'/t = 1, 
reaches a maximum value at an intermediate value of 
t'/t, and then drops to zero as t'/t — > 0. For instance, 
for n = 0.875, the maximum value AE S w 0.12t, which 
occurs for t'/t — 0.6, is approximately 4 times larger 
than its value in the uniform ladder. More broadly, we 
have studied the spin gap as a function of both U/t and 
t'/t; the results for n — 0.875 are shown in Fig. [3] One 
can see that AE S exhibits a broad maximum for U of 
order the band-width (U ~ 4 — St) and intermediate in- 
homogeneity, t'/t ~ 0.5. This figure looks qualitatively 
similar to the analogous result for the two dimensional 
checkerboard Hubbard model obtained previously by ex- 
act diagonalization of a 16 site system in Ref. [H; however, 
in contrast to that study, the present results are obtained 
in the thermodynamic limit. 

The dependence of AE P on U /t and t'/t is generally 
similar to that of AE S , as can be seen by comparing the 
contour plots of these two quantities for n = 0.875 which 
are shown in Fig. 2] and Fig. [3J respectively. However, 
there are interesting differences, as can be seen in FigJSJ 
where the two quantities are plotted as a function of t'/t 
for fixed U/t = 8 and n = 0.875. Note that for t'/t > 
0.6, AE p > AE S . This is, presumably, indicative of the 
existence of a spin 1 excitonic bound-state for t'/t > 0.6. 
A similar result was found in the uniform t— J mode l 13 ' 19 . 
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FIG. 2: The spin-gap, AE a , of the checkerboard-Hubbard 
ladder as a function of t'/t for n = 0.9375, 0.875, and 0.75 at 
fixed U = St. 
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FIG. 3: AE 3 of the checkerboard Hubbard ladder for n 
0.875 as a function of U and t', fixing t = 1. 



In order to calculate R = AEg / AE S , we must compute 
p c and £ s - To obtain £ s , we apply a relatively strong 
staggered Zeeman field of magnitude t to the end sites 
of the ladder and measure the decay of the staggered 
magnetization, M(j) as in Eq. [H] In all cases, we have 
found that M(j) decays rapidly on scales short compared 
to the length of our longest ladders, so £ s can be extracted 
from the calculations accurately. Representative results 
for M(j) are shown in the inset of Fig. [6] £, s as a function 
of t'/t is shown in Fig. [BJ for fixed n = 0.875 and U/t = 8. 
Note that for t'/t < 1/2, the spin-correlation length is 
roughly 3a, which is of the order of one unit cell of the 
checkerboard ladder. 

Next, we calculate both p c and K c following the pro- 
cedure described above [Eqs. [TOHT4"] . The value of K c is 
shown in Fig. for n = 0.75, 0.875, and 0.9375, fixing 
U/t = 8, as a function of t'/t. In contrast to the results 
for AE S (and somewhat to our surprise), for n = 0.875, 
K c is a weakly varying function of t'/t (and, as it turns 
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FIG. 4: The pair binding energy, AE P , of the checkerboard 
Hubbard ladder for n = 0.875 as a function of U and t , fixing 
t=l. 
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FIG. 5: AE a and AE P for the checkerboard ladder with fixed 
n = 0.875 and U/t = 8, as a function t'/t. 
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FIG. 6: The spin correlation length, £ s , for the checkerboard 
ladder with U = 8t and n = 0.875 as a function of t'/t, 
calculated from Eq. [9] The inset shows the expectation value 
of the spin (S z ) for U = St, n = 0.875, and t'/t = 1, as a 
function of position. A staggered Zeeman field of strength t 
has been applied to the two sites at the left edge of the ladder. 
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FIG. 7: Left: the Luttinger parameter K c as a function of 
t'/t for n = 0.75,0.875 and 0.9375 and U = St. The error 
bars were estimated by comparing between the values of K c 
obtained from Eq. [13] and Eq. 1141 Note that according to our 
definition of K c , the non-interacting value is K c — 2. Right: 
The phase stiffness p c (defined in Eq |12p as a function of t'/t. 



out, U/t as well). To a good approximation, for a wide 
range of values, we can simply take K c rs 1, independent 
of t'/t and U/t. Note that this implies that the super- 
conducting susceptibility diverges as T — > 0, so that it is 
reasonable to think of the ladder as a fluctuating super- 
conductor. (Of course, there is also a divergent charge- 
density wave susceptibility, xcdw ~ T~^ 2 ~ Ke \ so there 
is some unavoidable ambiguity with this simple intuitive 
picture.) As n is increased to 0.9375, K c increases, con- 
sistent with the expectation that K c — > 2 as n — > 

TABLE I: Values of the ratio R defined in Eq. [7]for n = 0.875 

and U = St. 

't' = 0.2 0.4 0.6 0.7 0.8 1.0 
R = 3.38 3.06 0.96 1.01 0.99 0.98 

From the measured values of £ s , k, and K c , the energy 
scale characteristic of phase-ordering can be extracted. 



Table U shows the ratio R from Eq. [7] Note that for 
t'/t > 0.5, R « 1. Thus, at least crudely, this regime 
can be thought of as a "BCS like" regime, in which there 
is a single energy scale, set by AE S , which characterizes 
the growth of superconducting correlations. Depending 
on precisely what criterion one chooses to quantify the 
crossover scale, phase fluctuations will produce a quanti- 
tative difference in the magnitude of the specified scale, 
but not large qualitative effects. Therefore, it is reason- 
able to assert that the values of the parameters which 
lead to the largest values of AE S and/or AE p are the 
"optimal values for superconductivity." 

For t' < 0.5t, we obtain R — 3, i.e. AEg > AE S , sug- 
gesting that this regime cannot be thought of in terms of 
either a naive weak or strong coupling picture. Remark- 
ably, the transition from R ~ 1 to R > 3 occurs quite 
sharply around t = 0.5, close to the point where the spin 
gap is optimal. 



The period 3 ladder 




t' (in units of t) 

FIG. 8: AE S for the three types of inhomogeneous ladders in 
Fig. [T]is shown at fixed n = 0.875 and U/t = 8 as a function of 
t'/t. The inhomogeneity induced by breaking up the ladders 
to period 1, 2 and 3 clusters increases the spin gap for t'/t < 1. 
The increase is most dramatic for the checkerboard ladder, in 
which the maximum spin gap is about 4 times larger than the 
spin gap for the uniform (t 1 — 1) system. For the period 1 
(dimer) ladder, the enhancement is by a factor of 2, while for 
the period 3 ladder the spin gap is only slightly enhanced (by 
about 10%). 

It is interesting to note that for n — 0.75, t'/t — 0.4 we 
find a sharp decrease of K c and p c . The value of K c at 
this point is smaller than the critical value of 1 at which a 
static charge-density wave should be stable^, indicating 
that this behavior of K c and p c may be due to a charge- 
density wave phase that exists for n = 0.75, t' ;$ 0.4i. 

We thus conclude that for the checkerboard Hubbard 
ladder, optimal superconductivity occurs for interme- 
diate values of U/t ~ 6, intermediate inhomogeneity, 
t'/t ~ 0.6 — 0.7, and electron concentrations near (but 
not equal to) one electron per site. We can now ask if 
this result is special to the checkerboard pattern, or if 
it applies more generally to the situation in which there 
are multiple sites per unit cell. We thus have repeated 
(although not in as much detail) the same calculations 
for the dimer ladder (period 1) and the period 3 ladder. 
(See Fig. [[J) 

In Fig. [8J we exhibit the dependence of the spin-gap 
of all three ladders for fixed U = 8t and n = 0.875 as a 
function of t'/t. In all the cases we see that there is an 
increase in the spin gap for some t'/t < 1. 

The result was expected, qualitatively, in the dimer 
(period one) case from previous works^^i, which found 
that the spin gap (as well as pairing correlations) is en- 
hanced upon making t' smaller than t in the dimer (pe- 
riod 2) ladder. In the case of the period three ladders, 
there is a very weak increase of the spin gap, which occurs 
at t'/t ~ 0.9. 

In Refs. [J, it was argued that the enhancement of su- 
perconducting correlations in the dimer ladder is due to 
the increase in the density of states close to the "Van 
Hove" point, in which one of the two bands of the two- 
leg ladder becomes unoccupied. Beyond this point, there 
is only a single band crossing the Fermi level, and the 
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FIG. 9: A system of coupled checkerboard ladders, connected 
with by a single particle tunneling matrix element t" 



system is likely to behave as a single-component Lut- 
tinger liquid. Therefore the superconducting signatures 
are rapidly suppressed. Consistently with this picture, in 
the dimer ladder, we find that the spin gap collapses to 
zero below t'/t < 0.6. In the period 2 (the checkerboard 
ladder) and period 3 cases, however, no such sudden sup- 
pression of the spin gap is observed as t'/t is reduced be- 
low the optimal point. This leads us to believe that the 
mechanism of the enhancement of the spin gap for t' < t 
in the checkerboard and period 3 ladders is unlikely to 
be related to a proximity to a Van Hove point. 

Note also that for all the inhomogeneous patterns in 
Fig. [51 the spin gap seems to reach zero at a critical 
t' c > (which is different for each pattern). In particular, 
for the "checkerboard" pattern, t' c ~ 0.05t. It is likely 
that for t' < t' c , the Luther-Emery phase gives way to a 
Luttinger liquid phase with one gapless charge mode and 
and gapless spin mode (or more) , although more work is 
needed to establish that. 

Overall, among all the patterns we have reported, the 
optimal ladder for superconductivity is a checkerboard 
ladder with U = 6t, t'/t = 0.6 - 0.8, and n = 0.875, for 
which AE S = 0.12i, AE p = 0.16t. 



VI. EXTENSION TO QUASI ID 

Above, we have argued that the superconducting ten- 
dency in the checkerboard-Hubbard ladder is optimized 
for an intermediate value of t'/t. However, since the 
superconducting T c of that system (as in any one- 
dimensional system) is strictly zero, one can worry that 
this statement may depend on how one chooses to mea- 
sure the strength of the superconducting correlations. 
We will now consider a system composed of an array of 
parallel checkerboard-Hubbard ladders coupled weakly in 
the direction transverse to the ladders, in which T c can 
be estimated in a controlled way based on the solution 
of the single-ladder problem. We will show that T c is 
maximal for ~ < 1. Thus, in this system, T c is indeed 
optimized when the electronic structure is non-uniform; 
i.e., there is an "optimal degree of inhomogeneity" for 
superconductivity. 

The quasi- ID system of coupled checkerboard Hubbard 
ladders is depicted in Fig. [S] The ladders are coupled by 
a single particle tunneling matrix element t" . We fix the 
value of t" <C t, t', and estimate T c (t'/t) from an inter- 
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chain mean field theory, described in Appendix A. From 
the numerical results for the checkerboard-Hubbard lad- 
der with n = 0.875 and U = 8t we recall that K c 

w 1 over the entire range < t' < 1 (see Fig. [7]). We 
therefore fix K c = 1, independent of t'. The resulting 
expression for T c is 

r ^'^W ^y. (15) 

Here, a; = v s /v c , and if (x) is a complete elliptic integral 
of the first kind. Note that T c depends on t'/t through 
v a , v c and AE S . As t' decreases, both v c and v s decrease; 
their ratio, however, is found to be approximately con- 
stant as a function of t'/t down to about t'/t = 0.5. 
(v s is obtained by using the estimate AE S £ S , where both 

AE S and £ s are calculated from DMRG.) AE S (y\ , on 

the other hand, has a maximum for ~ < 1. Therefore, 

as j is reduced from 1, T c (t') necessarily increases, and 
reaches a maximum for some t'„ < t. 



VII. DISCUSSION 

The present study, along with a variety of other re- 
cent studies£&2'22, provide strong support for a number 
of intuitively appealing ideas concerning the physics of 
the superconducting T c in unconventional superconduc- 
tors in which the pairing arises directly from the repulsive 
interactions between electrons: 1) The highest supercon- 
ducting transition temperatures occur at intermediate in- 
teraction strength, when U is comparable to the band- 
width. (A corollary of this is that materials which are 
studied because of their high transition temperatures are 
also likely to exhibit more general signatures of lying in 
an intermediate coupling regime; here, theoretical results 
from both weak and strong coupling approaches must be 
extrapolated, at best, to the limits of their regimes of 
applicability.) 2) Certain mesoscale structures ("optimal 
inhomogeneity"—) can lead to enhanced superconducting 
pairing, although clearly if the system is too strongly in- 
homogeneous, that always leads to a suppression of global 
phase coherence. 3) While short-range magnetic corre- 
lations, possibly of the sort envisioned in the putative 
RVB state of a quantum antiferromagnet or in certain 
theories of a spin-fluctuation exchange mechanism, may 
well be important for pairing, longer range magnetic cor- 
relations, especially of the sort one would expect near 
a magnetic quantum critical point, do not appear to be 
particularly favorable for superconductivity. (This final 
conclusion follows from a comparison of the t'/t depen- 
dence of the magnetic correlation length and the super- 
conducting pairing in Figs. |5]andEl respectively.) 

In addition, we found that the two-leg ladder at inter- 
mediate coupling (with U of the order of the bandwidth) 
and close to half filling is, in many respects, surprisingly 
well described as a "BCS-like" superconductor, in which 



there is a single crossover energy scale from the "nor- 
mal" to the "superconducting" state (rather than two 
separate scales, associated with pairing and phase coher- 
ence). This is based on the fact that the ratio of the 
pairing and phase coherence scales (defined in Eq. [7]) is 
close to its weak-coupling value, which justifies the iden- 
tification of the spin gap AE S (or the pair binding energy 
AE p ) as the relevant energy scale for superconductivity. 

Finally, there are a couple of unresolved issues and fur- 
ther directions we would like to highlight: 1) The extrap- 
olation of the present results to higher dimensions is, of 
course, the most important open issue. The strong qual- 
itative similarity between the present results and those 
obtained by exact diagonalization and CORE calcula- 
tions on relatively small 2D clusters certainly encourages 
us to believe that the results obtained here give insight 
into the behavior of the higher dimensional problem. In 
this context, it might be useful to carry out similar cal- 
culations on 4 leg and possibly even 6 leg ladders and 
cylinders, although it is considerably more difficult to 
extend these results to such long systems as are accessi- 
ble for the 2 leg ladder. 2) It is not clear exactly what 
aspects of the local electronic structure are essential fea- 
tures of an optimal inhomogeneity for superconductivity. 
In the present case, it is notable that pair-binding does 
not occur on an isolated dimer or six-site rectangle for 
any value of U/t, while there is pair-biding on an iso- 
lated square for U/t < 4.6. However, this observation 
does not provide an entirely satisfactory account of our 
findings, since the optimal pairing in the checkerboard 
ladder occurs for U/t = 4 — 8t, where the pair-binding 
energy of an isolated square is either small or negative. 
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VIII. APPENDIX: INTER-CHAIN 
MEAN-FIELD THEORY 

In this Appendix, we describe the inter-chain mean- 
field treatment of the quasi one dimensional system de- 
scribed in Sec. I VII This procedure is quite standard^— . 
We consider an array of plaquette ladders, modelled by 
Luther-Emery liquids. For simplicity, we will assume 
that each ladder is a single component system with a spin 
gap AE S . (The extension to the case of a two-component 
system is straightforward, and the result is qualitatively 
the same.) The ladders are coupled by an inter-chain 
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hopping term of the form: 

H± = -t" / dxip Pa (x, n) ip Pa (x,n + l), 

a,P=± n J 

(16) 

where i^ Pa (^j n ) (P — ±) creates a right or left mov- 
ing electron with spin a =f , i at position x in chain n. 



Next, we integrate out degrees of freedom of lengthscales 
smaller than the spin correlation length £ s ~ . Over 
such lengthscales, the system is essentially gapless and 
can be treated as a Luttinger liquid. To second order in 
t", the following effective inter-chain action is generated: 



Sf = (t") 2 J dxdr J dx'dT'{T^ +t<x (x, r, n) ^+, CT {x, t, n + 1) ^_ >tT , {x', r', n) ^_ iCr ; (a/, r', n + 1)) ,>, (17) 

aa' ,n 

i 



where (. . . )o,> denotes averaging over the "fast" (short- 
wavelength) degrees of freedom [using the decoupled 
(t" = 0) action], and T denotes time ordering. Since 
we are essentially performing a "coarse graining" step, 
increasing the cutoff of the theory from the lattice con- 
stant a to £ s , the region of integration in Eq. [17] is 

\J (x — x') 2 + v 2 (r — t') 2 < £ s . In order to evaluate 
the integrand, we write the fermionic fields in bosonized 



form: t(jp a ~ e i^{e a +Pip a ) ^ w h ere ^ an( j Q a are 
dual bosonic fields which satisfy [ip a (x) , d x a i {x')\ = 
id a(7 i5(x — x'). As usual, we introduce also charge and 
spin fields defined as ip c ,s = (ip^ ± ipi) /\/2, and similarly 
@cs = (^f±^4.)/v^' We define the fermionic Green's 
function Q (x, r) = (TV'+.f ( x j t j n ) $-,1 (0? 0; n ))o,> • Ex- 
pressing ^ (a;, r) in terms of the bosonic fields, 



G (x, t) ~ ( e 



V2T 2 ^ 2 1 2 1 2 



x 2 + (v c t) 



x 2 + (v c t) 



x 2 + {v s tY 



where we have used the shorthand notation 9 C = 9 C (x, r), 
6' c = 6 C (x', r'), 9 C = 9c (x/2, r/2), and similarly for 6 S , ip c 
and ip s . Plugging Q (x, r) into Eq. [TTland performing the 
integral, we get that the following inter-chain Josephson 
coupling term: 

Hf = -J ± J2 J dx$(x,n)& (x,n+l) + R.c, (19) 

where $ (x, n) = ipRfipn - ipRii>L^ ~ e iv/ ^ e <= cos V2ncp s 
and 




Here, X (a) = J^ 2 d\/\/l — a 2 sin 2 A is a complete el- 
liptic integral of the first kind. Eq. [20] contains a in- 
dependent pref actor, which we omit. 
The mean-field equation for T c is 

zJ ±X (Tc) = l, (21) 

where x (T) is the superconducting susceptibility of a sin- 
gle chain, and z is the number of nearest-neighbor chains 
(e.g., for a two dimensional array of checkerboard lad- 
ders, z — 2). Inserting Eqs. 12(201 in the mean-field equa- 
tion (|21l) . and using the fact that for the checkerboard 
Hubbard ladder K c w 1 over a wide range of parameters, 
we obtain Eq. [15] for T c . 
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